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ABSTRACT 



Aims. In this work we study the correlation between the soft (1.6-12.4 keV, mostly thermal) and the hard (20-40 and 60-80 keV, 
■ mostly non-thermal) X-ray emission in solar flares up to the most energetic events, spanning about 4 orders of magnitude in peak flux, 

$ZL(' establishing a general scaling law and extending it to the most intense stellar flaring events observed to date. 

Methods. We used the data from the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) spacecraft, a NASA Small 
Explorer launched in February 2002. RHESSI has good spectral resolution (^ 1 keV in the X-ray range) and broad energy coverage 
(3 keV-20 MeV), which makes it well suited to distinguish the thermal from non-thermal emission in solar flares. Our study is based 
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1 . Introduction erated to weakly relativistic energies and impact on the chromo- 

sphere, causing the heating and the evaporation of plasma into a 
Solar flares are violent explosions taking place in the solar coronal loop. The impulsive phase (lasting at most a few min- 
corona and chromosphere, generally attributed to magnetic re- utes in the Sun) is therefore characterized by hard X-rays of 
connection events. They result in the rapid release of a large non-thermal origins. The gradual phase (which in the Sun can 
amount of energy, up to 10^^ g^gs, in 10^-10^ s. last from a few minutes to several hours) is characterized by the 



^ ' on the detailed analysis of 45 flares ranging from the GOES C-class, to the strongest X-class events, using the peak photon fluxes in 

• the GOES 1 .6-12.4 keV and in two bands selected from RHESSI data, i.e. 20^0 keV and 60-80 keV. 

Results. We find a significant correlation between the soft and hard peak X-ray fluxes spanning the complete sample studied. The 
resulting scaling law has been extrapolated to the case of the most intense stellar flares observed, comparing it with the stellar 
^ ' observations. 

' Conclusions. Our results show that an extrapolation of the scaling law derived for solar flares to the most active stellar events is 

I compatible with the available observations of intense stellar flares in hard X-rays. 

I Key words, sun : flares - sun : X-rays - stars : flares - stars : activity - X rays - stars: individual : Algol, Ux Ari, AB-Dor 
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d A variety of phenomena is associated with the sudden energy soft X-ray emission which is (together with conduction to the 

■ ■ ■ release in a flare, including the acceleration of particles (elec- chromosphere) the main cooling mechanism of the evaporated 

Irons up to tens of MeV and ions up to tens of GeV) and the plasma filling a coronal loop, 
heating of coronal plasma to 10-30 MK. The radiation emitted 

by a solar flare covers the entire electromagnetic spectrum, from I" *is model only a small fraction of the energy of the accel- 

radio waves to X-rays and y-rays, but the diff-erent wavelengths crated electrons is lost through radiation: most of the loss is due 

are associated with diff-erent regions of the solar atmosphere and ^ Coulomb colhsions heating the chromospheric plasma, pro- 

with different mechanisms, making it necessary to take a "multi- coking its rapid heating at the loop foot point and leading to its 

wavelength" approach to the study of flare physics. evaporation (Antonucci et al. 1982, 1984; Fisher et al. 1985). 

In spite of the large number of solar flare observations per- The thick-target model implies a causal connection between 

formed with many instruments, the detailed physical model as- non-thermal (microwave and hard X-ray) and thermal (soft X- 

sociated with a flare remains elusive. The "non-thermal thick- ray and other) emissions, known as "Neupert effect." This is 

targetmodel" (Brown 1971; Hudson 1972; Lin & Hudson 1976), the correlation between the time-integrated non-thermal emis- 

in which a beam of electrons accelerated at the site of magnetic sion and the instantaneous thermal emission (Neupert 1968; 

reconnection (in the corona), are braked by the the cooler chro- Dennis & Zarro 1993; Veronig et al. 2002): as the instantaneous 

mospheric plasma (the "target"), producing the observed hard non-thermal radiation is a proxy for the amount of evaporating 

X-rays emission, has most commonly been used to explain many plasma, its time-integral should correlate with the total amount 

of the observed properties in solar flares. of thermal plasma, and thus radiation being observed in soft X- 

In this framework the flare has an impulsive phase and a rays, assuming that the time scale for cooling of the coronal 

gradual phase. In the impulsive phase the electrons are accel- plasma is much longer than the impulsive phase. 
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While the thick target paradigm has known difficulties by 
some authors, it has succeeded in explaining the emission ob- 
served in a number of very energetic stellar flares. In particular, 
the observed soft X-ray light curves of very large stellar events, 
lasting up to several days have been reproduced very well by 
numerical simulation of confined events in large coronal loops 
(Favata at al., 2005). The solar flaring mechanism seems then to 
be at work also in the most energetic stellar events, which can 
have peak energies of 4 to 5 orders of magnitude larger than so- 
lar events. 

In the present paper we address two questions: first, how 
does the hard X-ray component in solar flares (dominated by 
non-thermal radiation) coiTelate with the soft X-ray flux (domi- 
nated by thermal radiation), Dennis & Zarro 1993; Veronig et al. 
2002; Matsumoto et al. 2005. Second, if a correlation exists for 
solar flares could we also extend it to the case of the most inten- 
sive stellar flares? And by extension, can non-thermal emission 
be observed in stellar flares, with present, past or planned hard 
X-ray detectors? 

In contrast to a previous similar analysis by Battaglia et al. 
(2005), here we include the strongest X-class flares. A compari- 
son between the their and our results will be given in Sect. 5. 

The present paper is organized as follows: in Sect. 2 we give 
a brief description of the instrumentation used and the data anal- 
ysis; in Sect. 3 we present the results obtained for the case of 
solar flares; in Sect. 4 we extend these results to the case of stel- 
lar flares. Finally in Sect. 5 we discuss our results and in Sect. 6 
we summarize the main conclusions of the paper. 

2. Observations and data analysis 

RHESSI (Reuven Ramaty High-Energy Solar Spectroscopic 
Imager) is a NASA Small Explorer, launched on February 5, 
2002, which combines for the first time high-resolution imag- 
ing in hard X-rays and y-rays with high-resolution spectroscopy, 
in a way that a detailed energy spectrum can be obtained for the 
resolved solar disk (see Lin et al. (2002) and references therein 
for full details of RHESSI). 

The spectroscopy is achieved with nine cooled high-purity 
germanium crystals positioned behind the nine grid pairs of 
the telescope. These convert incoming X-rays and y-rays to 
charge pulses. The detectors have two electrically-independent 
segments: a front one to measure hard X-rays up to 200 keV and 
a rear one for energies up to 17 MeV, with spectral resolution of 
about 1 keV up to 100 keV and 3 keV up to 1 MeV. 

The RHESSI design features movable aluminum shutters to 
attenuate the fluxes of soft X-rays, preventing saturation and 
pulse pile-up that would occur at high count rates. The two shut- 
ter states used (Al, thin, and A3, both thick and thin in), are char- 
acterized by their half-response energies around 17 and 27 keV 
respectively. All of the analyses reported here are based on the 
simple diagonal elements of the spectral response matrix. While 
non-diagonal elements of the response matrix become important 
at the lower energies (i.e. bbelow 10 keV), we have only used, 
from RHESSI, data at £ > 20 keV, which justifies the approach 
chosen. 

2.1. Data selection 

Solar flares are classified on the basis of the GOES 
(Geostationary Operational Environmental Satellite) observa- 
tions (e.g., Garcia 1994), which monitor solar X-ray emission 
without spatial resolution (together with solar particle fluxes. 



etc.). The flare classification is based on the logarithm of the flux 
in the 1 .6-12.4 keV band, so that each flare is identified by a let- 
ter (as per Table[T]i plus a numeric code indicating the subclass. 



Table 1. The coiTespondence between the GOES class of a solar 
flare and the flux at Earth in the 1.6-12.4 keV band. 



GOES Class 


Flux 1.6-12.4 keV,W-m-^ 


A 


1.0 X 10-"* 


B 


1.0 X 10-'' 


C 


l.Ox 10-« 


M 


1.0 X 10^^ 


X 


1.0 X IQ-" 



Our aim was to study whether a general correlation exists 
between the peak soft and hard X-ray fluxes in solar flares, and to 
verify whether such coiTelation can also be extended to the case 
of intense stellar flares. We have therefore studied solar events 
spanning as broad a range of intensities as possible, from the 
weakest flares in the CI class up to the strongest available X 
class flares. 

We fixed a total number of around 50 events to study, and we 
tried to sample each class (C, M and X) evenly, with 15 events 
per class. The criteria used to select the flares in our sample 
were: 

- The flare must have good GOES data 

- The peak hard X-ray emission must be well observed (no 
RHESSI night or South Atlantic Anomaly passage near peak 
time) 

- The flare must be well defined and not occur during the de- 
cay phase of a larger one 

The list of flares which compose our sample is found in 
Table |2l Note that the quiescent emission has been subtracted 
from the values of the peak photon fluxes for all three bands 
of interest. The value of the quiescent emission has been deter- 
mined by selecting an interval in the pre-flare phase, possibly 
nigh time, for the RHESSI data and interpolated between the 
pre-flare and post-flare values for the GOES data. 

Many good candidates were available for class C or M is 
events, as thousands of flares have been observed in these two 
classes. The class C ad M events listed in Table |2] are therefore 
a representative sample, albeit by necessity somewhat arbitrary. 
On the other hand for the class X events the choice is limited, as 
the total number of class X events observed by RHESSI through 
2005 is only 62. Once our criteria are applied, 25 class X events 
are left, from which we selected our sample of 15 events. The 
clear correlation between the peak soft and hard X-ray emission 
found (e.g. Fig. [3]) confirms a posteriori that our events are rep- 
resentative and selected from a homogeneous parent sample. 

2.2. Data analysis 

For each flare in the list we produced light curves in a number of 
energy bands from the RHESSI data, using the RHESSI software 
(Schwartz et al., 2002) with a time bin of one rotation period of 
the RHESSI instrument collimator (^ 4 seconds). From the light 
curves in the relevant energy bands we determined the photon 
fluxes at the peak of the flare in the 20^0 keV and 60-80 keV 
bands. We performed a similar operation for the peak energy flux 
in the 1.6-12.4 keV band from GOES data. 
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Table 2. The list of flares used for our analysis. The time is the 
peak time in the 12-35 keV band as provided by the RHESSI 
Experimental Data Center (HEDC; Saint-Hilaire et al. 2002). 
The asterisks denote flares for which we were able to deter- 
mine the peak emission in both the 20^0 keV and in the 60-80 
keV bands. The classification under "GOES orig." is the original 
GOES class before subtraction of the quiescent emission, while 
the one under GOES is the one with the quiescent emission sub- 
tracted, which we used in our work. 



Flare Date Time GOES GOES orig. 



While the RHESSI passband in principle extends down to 3 
keV, we prefer to use the GOES data for the soft (mostly ther- 
mal) emission as the analysis of the RHESSI data in these bands 
is complicated by the presence of the moving attenuators for the 
more energetic events, as well as by the more complicated re- 
sponse of the detectors. 

We have chosen to compare the fluxes at the peak of the 
events, rather than e.g. total energy released in each band, as this 
is a quantity which can be determined more reliably in stellar 



flares, also for events with a modest signal-to-noise, for which 
one cannot reliably follow the flare decay (during which signifi- 
cant amounts of energy are still released). 

To make our analysis as non-parametric as possible, we have 
decided not to rely on spectral fits to the RHESSI data. To de- 
termine photon fluxes from the raw counts measured by the 
RHESSI detectors we have used the "semi-calibrated" data pro- 
duced by the RHESSI data pipeline. As described above this im- 
plies using of the diagonal elements from the response matrices 
of the detector For this analysis to be reliable, the following con- 
ditions have to be satisfied: 

1. The flare counting rate over the background should be at 
least ten times the background rate for the energies consid- 
ered. 

2. Only energies substantially below 100 keV can be dealt with 
properly in the semi-calibrated analysis. This is due of the 
increasing importance of Compton scattering at higher ener- 
gies. This has driven the choice of adopting 80 keV as the 
upper limit on our high-energy band. 

3. The fluxes below 10 keV are reliable only if no attenuators 
were in place. In the presence of attenuators the 3-10 keV 
band can be dominated by germanium K-shell escape events. 
This limitation in the RHESSI data has driven our choice to 
use data from GOES for the 1.5-12.4 keV band. 

We used only data from the front segments of the RHESSI 
detectors, as we focused on energies below 100 keV. We also 
discarded data from detectors 2 and 7 because of their bad reso- 
lution. 

The RHESSI pipeline software automatically coiTects for 
front-segment "decimation" (a RHESSI telemetry-saving fea- 
ture) which produced loss of counts below the decimation upper 
threshold energies. It also corrects for the presence of the attenu- 
ators. The pile-up coiTection is disabled by default but we enable 
this for all strong flares. 

Figs. [T] and 12] show the light curves of 2 flares in class C and 
M respectively. In addition to the light curves in the 20^0 and 
60-80 keV bands we also plot the emission in 3.-12.4 keV band 
from RHESSI, which for the weaker flares is not affected by the 
presence of attenuators and thus it is close to the GOES band 
1.6-12.4. We have however not used this in our analysis (relying 
on the GOES data). Note how the emission in the 60-80 keV 
band is barely visible (if at all) for the C class flare, while it 
becomes clearly present in the M class event. In both events the 
peak emission in the 20^0 keV band precedes the one in the soft 
component as expected if the hard emission is associated with 
the impulsive, non-thermal phase. The M-class event, by com- 
parison of the 60-80 and 20-40 keV RHESSI bands, shows the 
so-called "soft-hard-soft" spectral evolution as expected (Parks 
& Winckler, 1969; Hudson & Farnik, 2002; Fletcher & Hudson, 
2002; Grigis & Benz, 2005). 

3. Scaling laws for solar flares 

Fig. |3] shows the correlation between the peak emission in the 
1.6-12.4 keV band, as observed by GOES, and the RHESSI 20- 
40 keV band peak emission, for all flares in our sample (Table|2]i. 
The representative error bar on the RHESSI data has been esti- 
mated from the data itself combining the eiTor on the peak flux 
and on the background. For the peak we considered the fluctu- 
ations around the peak in a short binning of 4 seconds. Typical 
values were 15% for the 20-40 keV band and 10% for the 60-80 
keV band. For the background we took the peak to peak excur- 
sion and for the GOES band we estimated an eiTor of 20%. 
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The GOES flux vs the 20-40 keV HESSI flux 
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Fig. 3. The peak emission in the 1.6-12.4 keV band from GOES 
data as a function of the peak emission in the 20-40 keV band, 
from RHESSI data, for all flaring events in our sample. Blue 
Fig. 1. The light curves for flare 2 from Tabled which has GOES diamonds, green stars and red triangles correspond to C, M and 
class C2.7. Fluxes are in units of ph/(cm s keV). X class respectively. A representative error bar is shown in the 

lower right part of the plot. 



RHESSI Photon Flux vs Time 
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on the comparison between this scaling law and the one found 
in Battaglia et al. (2005) is included in the Section 5. 

In Fig.|4]we present the same type of analysis as in Fig.[3]but 
considering peak emission in the 60-80 keV band. Note that the 
number of events here is less than for the 20-40 keV band, as for 
the weaker flares no emission is detected in the 60-80 keV band. 
The flares for which a 60-80 keV peak flux has been determined 
are indicated with an asterisk in Table|2] 



The GOES flux vs the 60-80 keV HESSI 'lux 



Fig. 2. The light curves for flare 23 from Table |2] which has 
GOES class M5.6. Fluxes are in units of ph/(cm^ s keV). 
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The two quantities are well correlated, and a scaling law in 
the form of a power law can be derived, applying the bisector 
regression method as described in Isobe et al. (1990): 



0.10 1.00 
60-80 keV HessI 



10.00 100.00 1000.00 
flux in photons/(cm' s keV) 



0000.0' 



Fg ~ 7.83 X 10 



-6 pO.73 
^ (20-40) 



(1) 



where Fc is the GOES peak flux 1.6-12.4 keV band in units of 
Watt/m^ and F20-40 is the peak flux in the 20-40 keV band, in 
ph/(cm^ s keV). The slope for the power law is 0.73 + 0.04 and 
the intercept is (7.83 + 0.98) x 10"^. This correlation holds over 
more than 3 orders of magnitude in GOES peak flux. An analysis 



Fig. 4. The peak emission in the 1.6-12.4 keV band from GOES 
data as a function of the peak emission in the 60-80 keV band, 
from RHESSI data, for aU flare events wifli detectable 60-80 
keV fluxes. Green stars and red triangles correspond to M and X 
class events respectively. A representative eiTor bar is shown in 
the lower right part of the plot. 
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Although the spread is more significant than in the case of 
the 20-40 keV flux, a correlation between soft and hard peak 
fluxes is still well visible. We derived, using the same approach 
as for Eq.[ri a power law 

1.02 X 10-4 (2) 

where the units are the same in Eq.[T]and the slope is 0.58 + 0.07 
and the intercept (1.01 + 0.13) x 10"^. 



The relationship between the flux in the two bands is obvi- 
ously linear, with a proportionality constant which depends on 
the temperature. Writing the relation as 

^^(20-40) ^m-Fc (3) 

and using the same units as before, i.e. Fq in W/m' and f (20-40) 
in ph/ (cm^ s ke V), the values of m for different temperatures are 
given in Tabled for both the 20-40 and the 60-80 keV band. 



3. 1 . Thermal contribution to fhe hard band flux 

We have considered the total peak flux in each band, regardless 
of its origin. To determine which fraction of the observed flux 
is of thermal origin, we have developed a procedure to subtract 
from the total peak flux observed in the 20-40 and in the 60-80 
bands the thermal component, as a function of the peak flare tem- 
perature, using a set of model spectra. This is necessary to com- 
pare the predictions of our scaling laws with the observations 
of intense stellar flares, where the observed plasma temperatures 
are much higher than for solar flares and thermal emission can 
be significant also in these bands. 

To do this, we have simulated a set of thermal spectra and 
determined the relationship between the flux in the GOES 1.6- 
12.4 keV band and the flux in the 20-40 and 60-80 keV bands 
as a function of the peak temperature of the flare. We used 
flie XSPEC package (Schwartz 1996; Smith et al. 2002) with 
a MEKAL thermal model, a plasma emission code which imple- 
ments the optically-thin collisional ionization equilibrium emis- 
sivity model described in Mewe et al. (1995). For diff'erent tem- 
peratures we determined the value of the fluxes in the different 
bands of interest, 1.6-12.4, 20-40 keV and 60-80 keV. In Fig.|5] 
we plot the flux (as given by the model) in the 1.6-12.4 keV 
versus the flux in the 20-40 keV band for different temperatures 
ranging from 1 keV to 20 keV. 



The thermal 1.55 - 12.4 keV flux vs the thermal 20,-40.keV 
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Table 3. The values of m in giving the proportionality be- 
tween the thermal flux in the 20-40 and 60-80 keV bands 
(F(2o-40) and F(a)-s.O) and the thermal flux in the GOES 1.6-12.4 
keV band (Fc), with Fq in W/m^ and F(2o-4o) and F(6o-80) in 
ph/(cm^ s keV). 



Temp. (keV) 


m, 20-40 keV 


m, 60-80 keV 


2 


7.93 X 10* 


4.10 X lO-' 


3 


1.62 X 10* 


0.69 


6 


3.62 X lO'' 


1.27 X 10-* 


9 


1.08 X lO** 


3.63 X lO-'^ 


12 


1.92 X lO** 


1.98 X 10* 



While we cannot attribute a peak temperature individually 
to each flare in our sample, Feldman et al. (1995) found a tight 
correlation between peak flare intensity (up to GOES class XI) 
and peak temperature, allowing us to attribute an average tem- 
perature for some reference classes of flares. For each of these 
temperatures we calculate the thermal contribution with XSPEC 
and the corresponding total predicted flux obtained reversing the 
scaling law in EqlUwhich gives the following expression: 

^(20-40) ~ 9.908 X lOV^" (4) 

in the same units used above. 

We computed the ratio between the peak total flux (obtained 
from Eq. HJi and the average peak thermal flux in the 20-40 keV 
band. The results are shown in Table |4] 

Table 4. The ratio between the peak thermal and non-thermal 
flux in the 20-40 keV band as a function of the flare GOES 
class. The first column reports the class of flares, the second 
column the average peak temperature as derived from Feldman 
et al. (1995), the third column the peak thermal flux (in units 
ph/cm^s keV) from the MEKAL models, and the last column 
the ratio between the total peak flux derived by the scafing law 
in Eq.|4]and the average peak thermal flux. 



Class 


Average T (keV) 


.ftli(20-40) 


.f(20-40)/^'th( 20-40) 


CI 


1.30 


0.00058 


~ 100 


C4 


1.47 


0.02 


-24 


Ml 


1.64 


0.10 


~ 14 


M4 


1.70 


0.54 


~ 17 


XI 


2.00 


7.92 


^4 



Fig.5. The flux in the 1.6-12.4 keV band as a function of the 
flux in the 20-40 keV band for a MEKAL thermal models, for 
temperatures of 1, 2, 3, 6, 9, 12, 16, and 20 keV from left to right, 
which correspond respectively to ~12, 23, 35, 70, 104, 139, 186 
and 232 x lO^K. 



We can see that the thermal contribution to the peak flux in 
the 20-40 keV band is negligible in the less intense flares, and 
can be relevant in the 20-40 keV band only for the more intense 
events and and for temperatures close to the maximum that is 
observed on the Sun. 

Using the data in Table |3] and Eq. [T] for a temperature of 
around 2 keV, for instance, we find that a value above Fq ~ 3 x 
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10"^ W/m^ is required in order to have a total predicted F(2o-40) 
emission higher or equal to the thermal component. This means 
that temperatures of the order of 20 x lO^'K could be reached 
most likely for flares in M and X class, which is consistent with 
the result found by Feldman et al. (1995). 

In Fig. |6] we show the same data as in Fig. [3] for the peak 
emission during our sample of solar flares, with superimposed 
the relation between the thermal flux in the GOES band versus 
the thermal flux in the 20-40 keV band for a peak temperature 
of 1.5, 2.0 and 3.0 keV, spanning the range of temperatures ob- 
served in solar flares from the CI to the XIO classes. The class 
X flares are all located between the loci for the thermal fluxes 
for plasma at 2.0 and 3.0 keV (peak temperatures which are not 
exceptional for class X events). If the temperatures associated 
to our observed flares were lower than 2 keV, the fluxes would 
be totally dominated by non-thermal effects. For temperatures 
around 2.5 ke V a part of the emission in this band, for the intense 
solar flares, should be given by the tail of the thermal spectrum. 



The GOES flux vs the 20-40 keV HESSI flux 



The GOES flux vs the 60-80 keV HESSI flux 



10" 




0.10 1.00 10.00 100.00 1000.00 10000.0' 
20-40 keV Hessi flux in photons/(cm' s keV) 

Fig. 6. The same data as in Fig[3]with three additional dashed 
lines showing the thermal flux in the GOES band as a function 
of the thermal flux in the 20-40 keV band for peak flare temper- 
atures of 1.5, 2.0 and 3.0 keV (spanning the range observed for 
solar flares from the CI to the XIO class). 



In Fig. |7]we show the same data as in Fig. |6]for the 60-80 keV 
band. In this band the observed emission is entirely non-thermal: 
to contribute to the observed flux the thermal plasma would need 
to have temperatures ^ 6 keV, that is much higher than what is 
observed in solar events. 



4. Stellar flares 

In the previous section we have established a scaling law be- 
tween the peak flare emission in the GOES, 20-40 keV and 40- 
60 keV bands in solar flares, including the strongest X-class. In 
this section we go further and determine whether the same law 
extends to the emission from intense stellar flares. This is a step 
beyond all previous studies which provides interesting new re- 
sults. 
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Fig. 7. The same as Fig |4] with three additional dashed lines 
showing the thermal flux in the GOES band as a function of the 
thermal flux in the 60-80 keV band for plasma at 3.0, 6.0, 9.0 
keV. 

For this purpose we have compared the extrapolation of the 
solar scaling law to the fluxes observed in stellar events with the 
few available broad-band data from intense stellar flares. The 
only observatory which has to date performed broad-band ob- 
servations of stellar flares, from 0.1 up to ^ 100 keV has been 
BeppoSAX (Boella et al. 1997a). We have analyzed the data 
from the MFCS (1.6-10 keV) (Boefla et al. 1997b) and PDS 
(15-300 keV) (Frontera et al. 1997) instruments for the flares 
observed on UX Ari, (Franciosini et al. 2001), Algol, (Favata et 
al. 2001) and AB Dor (2 events, Maggio et al. 2000). For the 
purpose of the present paper we needed to derive the peak fluxes 
in the 3 bands of interest, distinguishing the total, thermal and 
non-thermal fluxes. 

We extrapolated the solar scaling laws, and used the peak 
flux in the GOES band (determined from our analysis of the 
BeppoSAX data) for each stellar flare to determine the predicted 
flux in the 20-40 and 60-80 keV band, using Eq.|4]and the fol- 
lowing one: 

/^(60-80)~ 7.680 xlO'^F^-^l (5) 

The results are shown in Fig. [8] for the 20-40 keV band and in 
Fig.|9]for the 60-80 keV band (where the stellar fluxes have been 
rescaled to a 1 AU distance, for ease of comparison). The hori- 
zontal error bars for the stellar events represent the uncertainties 
in the 20-40 keV flux resulting from the uncertainty in the slope 
of the best-fit correlation for the solar flares. 



Table |5]reports the results of this analysis for the 20-40 keV 
and 60-80 keV bands, respectively, giving, for each flare, the 
best-fit peak flux as deduced by the scaling law (Eq.|4]and Eq.|5] 
scaled to the respective distance of each star), with the uncer- 
tainty induced by the uncertainty in the best-fit scaling law. The 
best-fit for the total peak flux as measured on the BeppoSAX 
spectra is also reported. 

We have performed a new analysis of the data, and we tried 
different possible models to fit the spectra. In particular we in- 
vestigated if the presence of a non-thermal power law compo- 
nent could have produced a better fit, especially for the PDS in- 
strument. The results obtained showed that the precision of the 
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The GOES flux vs the 20-40 keV flux 




20-40 keV flux in photons/(cm^ s keV) 



Fig. 8. The scaling law for the peak emission during solar flares, 
extrapolated to the regime of intense stellar flares, showing the 
relation between the peak flux in the GOES band versus the 
peak flux in the 20-40 keV band. The solar correlation (Eq. |4|i is 
shown, extrapolated to the peak GOES flux of stellar flares. The 
diamonds, stars and triangles represent the data from C, M and 
X class solar flares on the bottom left side. The black crosses 
show the corresponding location for the flares observed in the 
named stars in the top right side, using their peak flux in the 
GOES band, with the horizontal error bar representing the un- 
certainty induced by the uncertainty in the slope of the best fit 
correlation. The diamonds on the left side of the data point for 
each star show the relation between the non-thermal peak flux as 
derived from the BeppoSAX data. 



data at hand is not sufficient to draw any conclusion. A multi- 
temperature fit or a one-temperature plus a power-law fit are 
equally statistically acceptable. This means that we cannot ex- 
clude for sure the fact that this component has already been ob- 
served with BeppoSAX but not even demonstrate its presence. 
A recent study by Osten et. al. 2007 also investigates possible 
non-thermal emission in a stellar flare. 

From the upper part of Table |5] it is evident that the exten- 
sion of the correlation observed for solar flares between the peak 
emission in the GOES band and the emission in the 20-40 keV 
band fits rather well the data for the intense stellar flares ob- 
served with BeppoSAX. In all cases, the observed peak flux in 
the 20-40 keV band is within the range of values predicted from 
the extension of the solar relationship. Given that this entails an 
extrapolation over the solar GOES band peak flux by four orders 
of magnitude, it is remarkable that the simple solar scaling law 
succeeds in predicting the 20-40 keV band flux also for stellar 
flares, also given that their peak plasma temperatures are much 
higher 

The lower part of Table |5] shows on the other hand that the 
extension of the solar correlation for the peak flux in the 60-80 
keV band does not do a good job at predicting the peak flux in 
stellar events. In all cases, the predicted 60-80 keV peak flux is 
much higher than the observed value, by one to two orders of 
magnitude. 

The correlation derived for the peak 60-80 keV flux in the 
solar case has a much steeper slope than the one derived for the 
peak 20-40 keV flux. In particular, when extrapolating toward 
very intense flares, the use of the solar correlations would pre- 



Table 5. For each of the 4 intense stellar flares observed by 
BeppoSAX, the peak flux in the 20-40 keV band and 60-80 keV 
as predicted by the extrapolation of the solar scaling law (column 
2, Eq.|4]and Eq.|5]l is given, together with its uncertainty (column 
3). In column 4 the peak flux determined from the observations 
(using a fit including a power-law component) is given. Units are 
ph/cm^ s keV. 



Event 


-f (20-40) 


Range 


^(20-40) 




scaling law 




observed 


Ux Ari,;;, 


0.47 X lO-* 


[0.15,2.12] X 10-" 


0.20 X 10-" 


Algolf 


2.86 X 10-" 


[0.89, 13.3] X 10-" 


1.86 X 10-" 


AB Dor'r 


0.64 X 10-* 


[0.23,2.5] X 10-" 


0.92 X 10-" 


AB Dor, 2^ 


3.35 X 10-" 


[0.11, 14.6] X 10-" 


2.28 X 10-" 




-f(60-80) 


Range 


-f(60-80) 




scaling law 




observed 


UX Ari 


1.81 X 10-" 


[1.60 X 10-\ 0.0053] 


0.16 X lO-' 


Algol 


1.29 X 10-3 


[1.06 X 10-", 0.043] 


1.34 X 10-^ 


AB Dor, 1 


1.42 X 10-* 


[1.61 X 10-^ 0.0039] 


0.23 X 10-5 


AB Dor, 2 


1.14 X 10-3 


[1.07 X 10-", 0.031] 


2.42 X 10-5 


" Franciosini et al. (2001) 






* Favataetal. (1999) 






^ Maggio et al. (2000) 






Maggio et al. (2000) 







diet that around class XIOOOOO (the most intense stellar flares 
are of class ^Xl 000000) the flux in the 60-80 keV band would 
exceed the one in the 20-40 keV band, leading to an "inverted 
spectrum" type of situation, contrarily to all available observa- 
tional constraints. 



The GOES flux vs the 60-80 keV flux 
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Fig. 9. The same as Fig. [8]but for the 60-80 ke V band. 



The results in Table |5] are obtained considering a model with an 
additional power law, whereas in Table. |6] we give the values of 
the peak fluxes considering just the contribution from the ther- 
mal component. 

From this table we can immediately deduce that in the 20- 
40 keV band almost all emission is of thermal origin while in 
the 60-80 keV, as for the Sun, it is the non thermal contribution 
which dominates. 
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Table 6. For each of the 4 intense stellar flares observed by 
BeppoSAX, the thermal contribution to the peak flux in the 20- 
40 keV band and 60-80 keV. Units are ph/cm^ s keV. 



Event 


pill 

^ (20-40) 


' (60-80) 


Ux Ari 


0.24 X 10-'' 


0.64 X IQ-' 


Algol 


1.2 X lO-'' 


9.6 X 10-' 


AB Dor, 1 


0.35 X lO-"* 


0.43 X 10-' 


AB Dor, 2 


1.4 X 10-"* 


4.60 X 10-' 



As the extrapolation to the intense stellar flares in this band 
of energy seems to not be in agreement with the data at hand, 
probably different mechanisms play a dominant role that be- 
comes more evident for the very intense flares at the highest en- 
ergies and that can justify the discrepancy in the extrapolation of 
the solar scaling law to stellar flares in the 60-80 keV band. 

Beyond the "standard paradigm" for solar flares, which we 
have identified with the thick-target model of Brown (1971), 
other possibilities exist. In particular coronal hard X-ray sources 
of the type observed by Frost and Dennis (1971) and Hudson 
(1978) now seem to be almost routinely observed in major events 
by RHESSI (Krucker et al., 2007). 

5. Discussion 

The correlation law in solar flares represented by our Eq. [1] 
can be compared to the result found in Battaglia et al. 2005. 
They studied a sample of RHESSI flares, by performing spectral 
fits to separate the thermal emission from the non-thermal one 
(modeled as a power law), and restricting themselves to flares 
up to GOES class M5.9 (no X-class events were considered). 
They correlate the flux in the GOES band with the non-thermal 
monochromatic flux at 35 keV, find the scaling law 

Fg = 1.8 X 10"V°f (6) 

where Fq is the flux in the GOES band and F35 is the non- 
thermal flux at 35 keV as obtained through individual spectral 
fits, in the same units as in Eq. [T] 

Considering the different approach adopted here, the two 
scaling laws are quite consistent, especially regarding the slope, 
in agreement with the expectation that the emission in the GOES 
band will be dominated by thermal emission for the complete 
range of temperatures spanned by solar flares, while at 35 keV 
the emission observed will be mostly of non-thermal origin. 

On the other hand caution should be used in comparing di- 
rectly Eq.[T]and Eq.|6]where different physical quantities appear 
We used integrated fluxes in the 20-40 ke V and 60-80 keV bands 
over all the energy, whereas Battaglia et al. used a monochro- 
matic spectrum at 35 keV. This may justify a steeper slope in 
Eq.|6]than the one derived by us in Eq.[T] 

We would like to stress the fact that we do not make spectral 
fits, which is also a major difference with the work by Battaglia. 
We decided to proceed in a non-parametric way because of the 
difficulty to estimate the low-energy cut-off of the power law 
from the data. This can produce a significant uncertainty in the 
determination of the thermal component. Data itself cannot con- 
straint well the two components and then the derived parameters 
could be very uncertain. 

In Fig. [To] we represent the same as in Fig. |9] with an addi- 
tional fit obtained combining the solar data and the stellar obser- 
vations. 



The GOES flux vs the 60-80 keV flux 




10-*^ - 



IQ-^ 10° 10" 10'° 10'^ 

60-80 keV Flux in photons/(cm' s keV) 

Fig. 10. The same as in Fig. |9] with an additional fit obtained 
combining the solar RHESSI data and stellar flares observations. 

This produces a scaling law 

F(6o-80) ~ 6.20 X l&Fl,^' (!) 

which coincidentally has the same slope as in Eq.|4]for the 20- 
40 band. We remind that the scaling law for the 60-80 keV band 
has been deduced from a limited number of flares in Table |2] as 
the less intense events have a negligible emission at these high 
energies. 

If future additional points would confirm the relation in 
Eq. |7] this would imply an universal slope across more than six 
orders of magnitude in flare strength, which will be a quite re- 
markable result. 

6. Conclusions 

We have studied the correlation between the peak flux emitted 
in three different energy bands (1.6-12.4 keV, or "GOES" band, 
20-40 keV and 60-80 keV) for a number of solar flares span- 
ning a broad range of peak flux in the GOES band (or "GOES 
class"). In the solar case, the GOES-band flux is almost com- 
pletely dominated by the thermal emission of the flaring plasma, 
while the 60-80 keV flux is entirely due to non-thermal emis- 
sion. This is likely due to the accelerated electrons which in the 
thick-target model of flares cause the evaporation of the plasma 
which is responsible for the thermal emission. The 20-40 keV 
emission contains contributions from both the tail of the thermal 
spectrum and the non-thermal, power law spectrum, with the for- 
mer increasing with the GOES class, which in solar flares is well 
correlated with the peak temperature of the flare. 

The analyzed sample of solar flares spans 3 orders of magni- 
tude in GOES class (from CI to XI 0), and both the 20-40 keV 
and the 60-80 keV peak flux are well correlated with the GOES 
band peak flux. The good correlation for the 60-80 keV band 
(which in the solar case is purely non-thermal in origin) may 
be construed to imply support for a "thick-target" type of frame- 
work, in which the thermal emission is linked to the non-thermal 
one by a causal relationship. 

We have then studied the applicability of these scaling laws 
to the case of intense stellar flares. Our sample is small, limited 
by the availability of the data: BeppoSAX. has been the only X- 
ray observatory with the broad-band response needed for these 
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observations, and only 4 events have been detected in the hard Proceedings of the Euroconference and lAU Colloquium 188, Santorini, 

bands. In the future, additional stellar flares may be observed Greece, June 2002, ed. H. Sawaya-Lacoste, ESA SP-505. 

with Suzaku, and, later with the proposed Simbol-X hard X-ray 

telescope. 

The extrapolation of the scaling laws derived for solar flares 
to the case of intense stellar flares (an extrapolation of 5 orders 
of magnitude in GOES peak flux) results in a quite remarkable 
good agreement between the predicted and observed 20-40 keV 
flux, while it over-predicts by some two orders of magnitude the 
60-80 keV flux. 

Future observations, increasing the statistic of intense flares, 
should allow to clarify the origin of this discrepancy. If a power 
law like the one in Eq. |7] is found, the extrapolation would be 
possible for both bands and the two scaling laws would be de- 
fined by the same slope, which would be an even more remark- 
able result. 
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